Universal geometric entanglement close to quantum phase transitions 
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Under successive Renormalization Group transformations applied to a quantum state \^) of finite correlation 
length ^, there is typically a loss of entanglement after each iteration. How good it is then to replace |^) by 
a product state at every step of the process? In this paper we give a quantitative answer to this question by 
providing first analytical and general proofs that, for translationally invariant quantum systems in one spatial 
dimension, the global geometric entanglement per region of size L ^ S, diverges with the correlation length as 
(c/12) log (^/e) close to a quantum critical point with central charge c, where e is a cut-off at short distances. 
Moreover, the situation at criticality is also discussed and an upper bound on the critical global geometric 
entanglement is provided in terms of a logarithmic function of L. 
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Introduction and aims.- Quantum phase transitions at zero 
temperature play a key role in the occurrence of important 
collective phenomena in quantum many-body systems. In this 
respect, a considerable effort has been applied throughout the 
last few years towards a theory of entanglement in extended 
systems, such as quantum lattice systems and quantum field 
theories. In particular, there have been several attempts to 
generalize the c-theorem of Zamolodchikov, which implies 
that the entanglement properties of a one-dimensional quan- 
tum system are somehow lost along Renormalization Group 
(RG) trajectories llll. A deeper understanding of this theo- 
rem using tools from quantum information science has shown 
that, under consecutive RG transformations, a translationally 
invariant quantum system in one spatial dimension may suffer 
from a monotonic loss of its amount of bipartite entanglement 
along the flow [2, 3], which can be explained in terms of a set 
of majorization relations (the so-called fine-grained entan- 
glement loss along RG flows [Ql)- This has in part motivated 
the application of a number of renormalization group ideas to 
novel representations of quantum states ||5J]. 

However, and in spite of the above findings, for a generic 
quantum system it is not yet known how close can its quantum 
state be to a globally separable state after each RG transfor- 
mation. In other words, the behavior along RG flows of the 
global entanglement in extended quantum systems — that is, 
the multipartite correlations that are shared by many differ- 
ent parties — still remains unclear in many aspects. So far, 
studies of global entanglement in extended systems have only 
been done for very specific quantum lattice models (see e.g. 
fiHIH]), and solid analytic results of wide generality are still 
missing. 

This letter deals with the above situation by providing first 
explicit analytical derivations, in the case of quantum systems 
in one spatial dimension and invariant under translations, of 
the universal properties along RG flows of the global geo- 
metric entanglement per region of size L, which we call £ 
lHlgt]. As we shall see, this measure of entanglement allows to 
quantify the fidelity between the quantum state and the clos- 
est separable state of contiguous blocks of size L along the 
flow. More precisely, we establish that near criticality and for 
one-dimensional quantum systems of finite correlation length 
^, the global geometric entanglement gets saturated when in- 



creasing the size L according to 
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close to a quantum critical point with central charge c IK 
where e is a regularization parameter at short distances that 
coincides with the lattice spacing for lattice systems. The 
above relation implies a logarithmic divergence of the satu- 
ration value of £ with the correlation length ^ when criticality 
is approached, so that the quantum system experiences a loss 
of multipartite entanglement along RG flows that decrease £,. 
Furthermore, the situation at criticality is also discussed, for 
which we provide the upper bound 



f < ^ loK ( - 



L 



(2) 



implying that the average geometric entanglement per block 
can not grow faster than a logarithmic function in the size L, 
in agreement with previous numerical estimations for bosonic 
and fermionic lattice models jsl]. 

Global geometric entanglement- To introduce the measure 
of entanglement that we use throughout this paper we ini- 
tially consider a pure quantum state of N parties \'^) £ Ti = 

(8)ili ''^'*'' where T^^'l is the Hilbert space of party i. Our 
aim is to quantify the global multipartite entanglement of | ^I*) . 
Following 121, this can be achieved by considering the maxi- 
mum fidelity | Amaxl between the quantum state |^) and all the 
possible separable and normalized states |$) of the N parties. 
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which is related to the minimum distance between \^) and the 
set of separable states. In order to have a measure of entan- 
glement that is zero when \^) is separable we take the natural 
logarithm, 



E{^) = -\og{\K 



(4) 



Here we will be interested in the above quantity per party, 
which has a well defined infinite-A^ limit: 



fjv = N-^E{-^) , £ = lim £_ 



N- 



N 
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In the case of a one-dimensional quantum system of infinite 
size we choose the different parties to be contiguous blocks 
of size L as shown in Fig.([T]i. In this setting, £ corresponds 
to the global geometric entanglement per region of size L [8t|, 
and is the entanglement measure that consider in this work. 
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Figure 1: An infinite one-dimensional quantum system is divided 
into different parties corresponding to contiguous blocks of size L. 

Loss of global entanglement along RG flows.- Our aim now 
is to provide a derivation of Eq.([T]i in the introduction. The 
complete proof of that relation depends on a number of results 
that we have derived by using properties from linear algebra 
and conformal field theory. In the proof presented here, some 
strong technical details are avoided and referred correspond- 
ingly, always looking for the clarity of expression. 

To start with, let us consider the quantum state |^). For the 
moment, the total length of the system is assumed to be NL, 
and the limit N ^ oo will be taken shortly. The first step 
in our derivation is to find a suitable decomposition of state 
1 4') as in illH : first, we consider the bipartition of the system 
[1 : 2, . . . , TV] and compute the Schmidt decomposition of 
I*), 
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components M^ 
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In the limit N ^ oo and for a 
system in one dimension invariant under translations divided 
in blocks (parties) of equal size L, we make three natural as- 
sumptions about the behavior of the physical system: (i) it 
can be correctly described by a decomposition like the one in 
Eq.®, (m) its local description is site-independent, and (iii) 
the maximization from Eq.© can be done with a state |$) that 
is the tensor product of the same state \(f)) for all the parties. 

In this situation, we have that |0W) ~ |0), |t^U) = l^ap), 
A/[*l = M and /?['! = D for every party i. The fidelity |A| is 
then given by 



|A| = lim \d\ 



N 



(10) 



where d is the eigenvalue of largest magnitude of the matrix 
VDMVD. Notice that the fidelity |A| is zero for an infinite 
system unless it is a separable state of the blocks (see e.g. 
il2ll and references therein). However, the global multipartite 
entanglement per block of size L is finite and reads 



£ = -log(|dn 



(11) 



where |dmax| is the maximum possible eigenvalue of ma- 
trix ^/DM^/D over all possible quantum states 10) such that 
{(j)\(j)) = 1 (which appear in the definition of matrix M). For a 
given state |0), the maximum eigenvalue of matrix \fDM\fD 
can be obtained by solving a different maximization problem, 
namely. 



|dp = [max r t ^/dM^/d i 



(12) 



In the above equation, AL are the Schmidt coefficients, and 

I Tq ) , |tjQ ' ' ) are the left and right Schmidt vectors respec- 
tively. Next, we find the Schmidt decomposition of quantum 
state |tjQ ' ' ') according to the bipartition [2 : 3, 



,Nl 



/3 


'"') 


SO that state ^) reads 
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Proceeding iteratively as above for all the parties, the quantum 
state I '$) is finally expressed in terms of the Schmidt coeffi- 
cients for different contiguous bipartitions as IISJ I 



i*)= Y. AWA^••A 
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The above decomposition allows to obtain a useful expres- 
sion for the fidelity |A| ~ |($|^)| between |^) and some 
separable state 1$) = |0[il)|(?!)[2l) • • • l^'^') of the N parties, 
in the limit N ^ oo and for a translationally invariant state. 
To see this, let us previously define D^'' as the diagonal ma- 
trix of components D^^L = \aSaf:), and M^'l as the matrix of 



where f^f — 1. Therefore, to find |dmax| we need to solve 
two different maximization problems: one over the quantum 
states \(f>) and another over the vectors r. 

Our interest is now focused on solving this double max- 
imization problem. In order to achieve this, we fix vector 
r and perform the maximization over the quantum state 1 0) 
only. As a result, we obtain that the optimal state |(/)max) is 
given by |0,nax) == \'4'{r))/ Vi^-'i^Mi^)' with \^p{r)) = 

We then have that 




max |('(/;(r)|?/;(r))| 
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In the above expression we have introduced matrix A{L) for a 
block of size L. The components A{L)i^aa').{i3/3') of this ma- 
trix (where {aa') is understood as a single index and similarly 
for {(3(3')) are computed as shown in the diagram of Fig.©. 

Before moving to the maximization of Eq.(fT3]), we wish 
to explain certain spectral properties of matrix A{L) that 
turn out to be of relevance for our purposes. First of all, 
notice that A{kL) = A''{L) for k = 2,3..., see Fig.©. 
The reason for this is that a block of size 2L can be under- 
stood by grouping together two contiguous blocks of size L, 
which at the level of the matrix involves the multiplication 
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Figure 2: (color online) Diagrammatic representation of the compo- 
nents A(L)((j(j/) (^^/) of matrix A{L) for a block of size L. In the 
diagram, the violet diamonds correspond to the components ^fDafi 
of matrix \AD and the half-ellipses correspond to quantum states 
|ra/3). The scalar product between two quantum states is represented 
by two half-ellipses together, one of them representing the bra { \ 
and the other the ket \ ) . The different emergent lines correspond to 
the Greek indices a, l3..., where common legs between objects cor- 
respond to shared summed indices, and free legs conespond to free 
indices. 



of A{L) by itself. Moreover, we have the chain of equalities 
1 = (^1^) = limfc_ootr(A(L)^) = limfc^oo {i^i{A{L)))\ 
where i'i{A{L)) is the eigenvalue of largest magnitude of 
A{L). These relations necessarily imply that i'i{A{L)) = 1 
and |i^i(A(L))| < 1 for the rest of the non-zero eigenval- 
ues i'i{A{L)), i = 2,3... of A{L). In fact, there is only 
one eigenvalue of magnitude 1 in the case of systems away 
from criticality. This is so since the correlation length ^ of 
a one-dimensional quantum system in state |^) is given by 
(^/e) = — 1/ log |i^2(^(e))|, where e is a cut-off parameter at 
short distances that coincides with the lattic e sp acing in the 
case of systems defined on a lattice (see e.g. Ill3ll for a deriva- 
tion of this property). 




V^ 



A{L) 



A{2L) = A\L) 



_> 



A(3L) = A\L) 



Figure 3: (color online) The repeated multiplication of matrix A{L) 
for a block of size L produces the same matrix but for blocks of larger 
size, i.e. A{kL) = A^{L), for fc = 2, 3.... 

As a consequence of the above facts, a key property of ma- 
trix A{L) can be derived: the spectral decomposition of A{L) 
reads 



^(L)=5^(^,(A(6)))(^/^)a(z)a(z)t 



(14) 



where a{i) is the eigenvector corresponding to the ith eigen- 
value of largest magnitude. From this equation it is pos- 
sible to see that, if L ^ ^, then all the eigenvalues for 
i > \ are exponentially suppressed as exp{—{L/^i)) with 
(^,/e) = \og\vi{A{e))\, and therefore A (i) ~ a(l)a(l)t 
since only the eigenvalue of largest magnitude vi{A{e)) = 1 
does not vanish. Under these circumstances, the maximiza- 
tion from Eq.lfTSll transforms into 



max (r ( 



)Ml)| 



(15) 



which is the maximization of a scalar product between two 
vectors. 

In order to find |dmax| in Eq.dTSll. we make use of some 
properties of the decomposition in Eq.(|9]) of state |^). In 
particular, we use the fact that the corresponding eigenvec- 
tor a(l) has components a{\)aa' ~ ^a^aa', which is a con- 
sequence of the orthonormalization of the different left and 
right Schmidt vectors in Eq.(|9l) Ill9ll . The optimization from 
Eq.dTSll then gives IdmaxP = '^i, where Ai is the largest 
Schmidt coefficient in the decomposition of \'^) given in 
Eq.(|9l). The global geometric entanglement per block of size 
L is then given by 



f = -21og(Ai) 



(16) 



which holds away from criticality for L 3> ^. 

At this point the global entanglement £ only depends on 
the largest Schmidt coefficient Ai. However, it is possible to 
obtain a more convenient expression for £ in terms of the cor- 
relation length ^. This can be achieved by considering the 
reduced density matrix p of one of the blocks, given by 



^(AaA/3) \Tap){Tal3\ 
a, (3 



(17) 



(see e.g. 111411 for details on this derivation). Remarkably, the 
following inequality holds: 



i'yi{p)y/'>xi>iMp)+Mp)y^' 



(18) 



where J^i(p) refers to the ith largest eigenvalue of the re- 
duced density matrix p, and t^„(p) refers to the smallest 
one 1231 ■ Given the fast decay of the eigenvalues of p 
away from criticality (see e.g. the second reference in ^ 
for a discussion about this property), we can safely assume 
that Vn{p)/v2{p) ~ 0. Also, following the results from 
Sec. IV in yfl, it is not difficult to see that — log(i^2(p)) ^ 
— log {vi (p)) ^ (c/6) log (^/e) if the system is close enough 
to a quantum critical point with central charge c (so that ^ is 
large), which implies our claim in Eq.([T]) in the introduction. 
Upper bound for critical systems.- At a quantum critical 
point, the correlation length S, of the system diverges, so that 
several of the eigenvalues of matrix A{L) in Eq.(fT4l) are ex- 
pected to be of magnitude one and not necessarily real. This, 
in turn, makes the maximization of Eq.(fT3]) quite difficult. 
However, it is still possible to derive a general upper bound 
on the scaling of £ with the size L of the blocks. 



This upper bound is derived as follows: for a finite sys- 
tem of M parties, the squared fidelity |Ap between a mixed 
state pm for the M parties and a separable product state 
1$) = 10)®^ is given by |A|2 = ($|pjv|$). We can as- 
sume that the reduced density matrix pM describes the de- 
grees of freedom of M contiguous blocks in a pure quantum 
state of N > M parties, which is described by the decompo- 
sition in Eq.(|9]l. The density matrix pM is then obtained by 
tracing out the degrees of freedom of the remaining N — M 
sites in Eq.©. If the whole system is translationally invariant, 
and for infinite N, it is possible to define a site-independent 
vector b of components fe(QQ') = ^/\^Mal3^f\^ M^^, pX-^ 
(where matrix M is defined as in the off-critical case). The 
squared fidelity |Ap is then expressed in terms of vector b as 
|A|2 = fot {A{L))^'^'^b, where matrix A{L) is defined as 
previously. Next, we make use of the normalized vector a(l) 
of components a(l)(„Q/) = AqJqq'. Introducing a resolu- 
tion of the identity operator before and after each one of the 
matrices A{L) in the expression for |Ap, and in such a way 
that the resolution includes the projector a(l)a(l)t, we obtain 

|Ap = (6ta(^in +T, where T is some positive constant. 
It is possible to see that this expression is indeed equivalent to 



bounded in the limit i\f ^ cxj as 



|A|^ = ((0|p|</.))-+T, 



(19) 



where p is the reduced density matrix of a block of size L. 
From the above relation, we derive after some manipulation 
the inequality |A,„axP > {vi{p)Y^ for the maximum over- 
lap |A,„axP, where vi{p) is the largest eigenvalue of p. The 
global geometric entanglement £ per block is then seen to be 



£ < -\og{Mfi)) = E^{p) , 



(20) 



where Ei (p) is the single-copy entanglement between a block 
of size L and the rest of the system. Now we make use of the 
existing results on the behaviour of Ex (p) at the critical point 
of a quantum phase transition from {VJA- In particular, we use 
the property that, at criticality, the single-copy entanglement 
Ei{p) scales with the size L as Ei{p) — (c/6) log (L/e) for 
large L, where c is the central charge of the underlying con- 
formal field theory in (1 + 1) dimensions. By combining this 
result with Eq.(l20li. the expression from Eq.© in the intro- 
duction follows, providing a universal upper bound on £ for 
any critical quantum system in one spatial dimension in terms 
of a logarithmic function of the size L of the blocks. 

Conclusions.- Here we have established first analytical 
derivations of the global geometric entanglement per block 
of size L for quantum systems in one spatial dimension 
and invariant under translations. We have proven that one- 
dimensional quantum systems tend to be globally separable 
along RG flows by following a universal scaling law in the 
correlation length ^ of the system. Furthermore, an upper 
bound on the critical scaling of the global geometric entan- 
glement has been given in terms of the logarithm of the size 
L. Our results are analytical, universal, and are valid for all 
one-dimensional quantum systems close to and at criticality. 
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